Observational Constraints on Visser's Cosmological Model 
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Theories of gravity for wliicli gravitons can be treated as massive particles have presently 
been studied as realistic modifications of General Relativity, and can be tested with cosmological 
observations. In this work, we study the ability of a recently proposed theory with massive gravitons, 
the so-called Visser theory, to explain the measurements of luminosity distance from the Union2 
compilation, the most recent Type-la Supernovae (SNe la) dataset, adopting the current ratio of 
the total density of non-relativistic matter to the critical density (f2m) as a free parameter. We 
also combine the SNe la data with constraints from Baryon Acoustic Oscillations (BAO) and CMB 
measurements. We find that, for the allowed interval of values for Qrn, a model based on Visser's 
theory can produce an accelerated expansion period without any dark energy component, but the 
combined analysis (SNe la-f-BAO-l-CMB) shows that the model is disfavored when compared with 
ACDM model. 

PACS numbers: 98.80.-k, 95.36.-fx, 95.30.Sf 



I. INTRODUCTION 



The current Universe's energy budget is a consequence 
of the convergence of independent observational results 
that led to the following distribution of the energy 
densities of the Universe: 4% for baryonic matter, 23% 
for dark matter and 73% for dark energy [l|. The 
key observational results that support this picture are: 
mesurements of luminosity distance as a function of 
redshift for distant supernovae [2-4], anisotropics in the 
Cosmic Microwave Background (CMB) observed by the 
WMAP satellite Pi and the Large Scale Structure (LSS) 
matter power spectrum inferred from galaxy redshift 
surveys such as the Sloan Digital Sky Survey (SDSS) [Q] 
and 2dF Galaxy Redshift Survey (2dFGRS) Q. 

In order to explain all the currently available 
cosmological data, the cosmological concordance model 
ACDM need to appeal to two exotic components, the so 
called dark matter and dark energy. The latter drives the 
late time accelerated expansion of the Universe and it is 
one of the greatest challenges for the current cosmology. 
Indeed, the physical nature of the dark energy is a 
particularly complicated issue to address in the ACDM 
context, due to its unusual properties. It behaves 
as a negative-pressure ideal fluid smoothly distributed 
through space. One can ask if the accelerating expansion 
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of the Universe might indicate that Einstein's theory 
of gravity is incomplete, i.e., can an alternative theory 
of gravity explain consistently the late-time cosmic 
acceleration without recurring to dark energy? 

There are several alternative approaches based on the 
idea of modifying gravity. Currently, one of the most 
studied alternative gravity theories is the so called f{R) 
gravity, whose basic idea is to add terms which are powers 
of the Ricci scalar R to the Einstein- Hilbert Lagrangian 

Recently, M. Visser proposed a modification of the 
general relativity (GR) where the gravitons can be 
massive particles [IJ] . In particular, several authors have 
studied the limits that can be imposed to the graviton 
mass using different approaches. For example, from 
analysis of the planetary motions in the solar system 
it was found that rrig < 7.8 x 10~^^g [li|. Another 
bound comes from the studies of galaxy clusters, which 
gives rrig < 2 x 10~^^g [la). Although this second limit 
is more restrictive, it is considered less robust due to 
uncertainties in the content of the Universe in large 
scales. Studying rotation curves of galactic disks, de 
Araujo and Miranda [131 have found that rrig <^ 10~^^g 
in order to obtain a galactic disk with a scale length of 
6-10 kpc. 

Studying the mass of the graviton in the weak field 
regime Finn and Sutton have shown that the emission 
of gravitational radiation does not exclude a non null 
(although small) rest mass. They found the limit mg < 
1.4 X 10~^^g [i3l analyzing the data from the orbital 
decay of the binary pulsars PSR B1913+16 (Hulse- Taylor 
pulsar) and PSR B1534+12. 

In particular, as discussed by Bessada and Miranda 
[13 l. if rrig > 10~^^g then massive gravitons would leave 
a clear signature on the lower multipoles (/ < 30) in the 



cosmic microwave background (CMB) anisotropy power 
spectrum. Moreover, massive gravitons give rise to a non- 
trivial Sachs- Wolfe effect which leaves a vector signature 
of the quadrupolar form on the CMB polarization [20| • 

An interesting result that comes from Visser's model 
is that the gravitational waves can present up to 
six polarization modes |2l| instead of the two usual 
polarizations obtained from the GR. So, if in the future 
we would be able to identify the gravitational wave 
polarizations, we would impose limits on the graviton 
mass by this way. 

The Visser's theory of massive gravitons can be used 
to build realistic cosmological models that can be tested 
against available observational data. It has the advantage 
that it is not necessary to introduce new degrees of 
freedom neither extra cosmological parameters. In fact, 
the cosmology with massive gravitons based on the 
Visser's theory has the same number of parameters of 
the flat ACDM model but no extra fields are added. 
In this paper we derive cosmological constraints on the 
parameters of the Visser's model. We use the most recent 
compilation of Type-la Supernovae (SNe la) data, the 
so-called Union2 compilation of 557 SNe la [24]. We 
also combine the supernova data with constraints from 
Baryon Acoustic Oscillations (BAO) [H and CMB shift 
parameter measurements |24| . 

The paper is organized as follows: in Section |TT] we 
briefly review the Visser's approach. Section IIIII is 
devoted to the description of the cosmological model. In 
Section IIVI we investigate the observational constraints 
on the Visser's cosmological model from SNe la, BAO 
and CMB shift parameter data. In Section |V] we present 
our conclusions. 



II. THE FIELD EQUATIONS 

The full action considered by Visser is given by [14] 
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where besides the Einstein-Hilbert Lagrangian and the 
Lagrangian of the matter fields we have the bimctric 
Lagrangian 

Cmass{g,go) = lni'V^oUgo'r{g-go)^.AgoTi2) 



^{g- go) p. -I [{go Tig -go) 



fii^i 



where G'"^ is the Einstein tensor, T'"^ is the energy- 
momentum tensor for perfect fluid, and the contribution 
of the massive tensor to the field equations reads: 



AfM^ = (ffo-i)^^ 



(3-3oW-i(<?o).p(5o"^r^ (4) 



Note that if one takes the limit TOq 



(go) 



l\pu 



the usual 
Einstein field equations are recovered. 

Regarding the energy-momentum conservation we will 
follow the same approach of [25] and (2al in such a way 
that the conservation equation now reads [23, [23] : 



V^T^"" 



mTrGff 



V^M'^^ 



(5) 



since the Einstein tensor satisfies the Bianchi identities 



III. COSMOLOGY WITH MASSIVE 
GRAVITONS 

For convention we use the Robertson- Walker metric as 
the dynamical metric: 



ds^ = c^dt^ -a^it) 



dr^ 



1 — kr^ 



+ r^{de^+sm'^0d<p^) 



(6) 

where a{t) is the scale factor. The flat metric is written 
in spherical polar coordinates: 



dsl = c^df - [dr 



'■edcj,')]. (7) 



The choice of Minkowski as the non-dynamical 
background metric ga is based on the criterion of 
simplicity. In first place, the metric go is defined in 
such a way that it coincides with the dynamical metric 
g in the absence of gravitational sources. The other 
point is that we do not need additional parameters for 
the cosmological model. The last important point is 
that considering Minkowski for 170 we obtain a consistent 
relation for the energy- momentum conservation law [27| . 

Using © and ([7]) in the field equations ^ wc get the 
following equations describing the dynamics of the scale 
factor (taking A; = for simplicity): 



1 



-TO^c^(a^ - 1) = 



SttG 



(8) 



where m = jrigc/h, nig is the graviton mass and {go)tiii 
is a general flat metric. 

The field equations, which are obtained by variation of 
([1]), can be written as: 



2 c4 
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and 



-m c a [a — 1) = — 
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where as usual p is the energy density and p is the 
pressure. 



From Eq. ([5]) we get the evolution equation for the 
cosmological fluid, namely: 
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0, (10) 



where H — a/a. Considering a matter dominated 
universe (j) = 0) the above equation gives the following 
evolution for the energy density: 
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where po is the present value of the energy density. Note 
that in the case rUg — j' we obtain the usual Friedmann 
equations. 

Now, inserting dill) in the modified Friedmann 
equation ([8]) we obtain the Hubble parameter: 



H\a) = Hi 



nl 



1 



+ ;-nl {7a' - 5a') 



(12) 



where the relative energy density of the z-component is 
^i = Pi/Pc {pc = 3H'c'/8ttG is the critical density) 
where '«' applies for baryonic and dark matter. Moreover, 
the present contribution of the massive term is defined 
by: 



n^„ 



1 

70 






(13) 



where mn = TiH^^/c' is a constant with units of mass. 

Since we are assuming a plane Universe (fc = 0), the 
total density parameter is ^lotai — 1- Thus, V>P can be 
replaced by 57^ = 1 — J7°n- This tell us that the model 
described by the Hubble parameter ([T^ has only two free 
parameters, namely Hq and fij^j, which can be adjusted 
by the cosmological observations, i.e., the same number 
of free parameters of the ACDM model. 



Also, fJ-obs{zi) are the values of the observed distance 
modulus obtained from the data and a{zi) is the 
uncertainty for each of the determined magnitudes from 
supernova data. 

Evaluating the minimum value of x^ from the Union2 
compilation of SNe la [22>| we found Xmin — 561.11 for 



the Visser's theory, with il^. 



0.261] 



have considered errors at 1 sigma level. 



0.020' 



where we 



B. Baryon Acoustic Oscilations 

The primordial baryon-photon acoustic oscillations 
leave a signature in the correlation function of luminous 
red-galaxies as observed by Eisenstein et al. i23']. This 
signature provides us with a standard ruler which can be 
used to constrain the following quantity 
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dz 
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where E{z) = H{z)/Hq, the observed value of A is 
Aobs = 0.469 ±0.017 and zi = 0.35 is the typical redshift 
of the SDSS sample. The computation of the value of ftm 
which better adjust Aobs lead us to f7,„ = 0-306]']g'Q25- 



C. CMB Shift Parameter 

The shift parameter R, which relates the angular 
diameter distance to the last scattering surface with the 
angular scale of the first acoustic peak in the CMB power 
spectrum, is given by (for fc = 0) [2J, 1291 



1089 



R 



1089 



n, 



772 



dz 

^ = 1.70 ±0.03. (18) 

H 



IV. ANALYSIS AND DISCUSSION 



A. Supernova la 

In order to put constraints on the cosmological model 
derived from the Visser's approach, we minimize the x^ 
function 

^2^^^^^y [Pt,{z^\n^)~Pobs{z.)f 
^-^ (y^\Zi) 

where /it/i(zi |r2,„) is the predicted distance modulus for 
a supernova at redshift zi. For a given Vim we have 

H{z\n,n) = m~ M ^ 25 + 51ogdL(z|ri,„), (15) 

where m and M are, respectively, the apparent and 
absolute magnitudes, and dL{z\^m) stands for the 
luminosity distance given by 



dL(z|a«) = (l + z)c 



dz' 



H{z'\n^y 



(16) 



It is worth stressing that the measured value of -Rioso is 
model independent. Also, note that in order to include 
the CMB shift parameter into the analysis, it is needed 
to integrate up to the matter-radiation decoupling [z ~ 
1089), so that radiation is no longer negligible and it was 
properly taken into account. With these considerations, 
the best-fit value for the relative matter density using 
_Ri089 is i^m = 0.224_g;Q3g. 



D. Joint analysis 

When the measurements of SNe la luminosity 
distances are combined with information related to the 
Baryon Acoustic Oscillation (BAO) peak and the CMB 
shift parameter, the constraining power of the fit to 
the parameters in the cosmological model is greatly 
improved. Following such an approach we examine here 
the effects of summing up the contributions of these last 
two parameters into the x^ of Eq- (HH)- Our result is 
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TABLE I: Best-fit values for Q.m for tlie cosmological observables considered in tliis work. It is also shown how the introduction 
of systematic errors from the SNe measurements can affect the best-fit. We have worked only with fiat Universe models, i.e., 
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FIG. 1: Hubble parameter as a function of the redshift for 
best-fit value obtained from SNe la. By using the different 
best-fit values, the curve does not change significatively 



^m = 0.273 ± 0.015 with the corresponding minimum 
value for the y^ function: Xmin — 565.06. 

We can compare our resuhs with the ACDM model by 
taking the difference between x^q ^-nd Xacda/' "which are 
the minimum x^ values for the massive bimetric model 
and for the ACDM model, respectively. The evaluation 
of this difference gives the result Ax^ ~ ^g ~ Xacdm = 
21.30, which shows that the bimetric Visser's model is 
disfavored when compared with the flat ACDM model. 

In the Table U we summarize our results for 51^ 
considering each cosmological observable: SNe, CMB, 
BAO and the combined analysis (SNe-^CMB-^BAO). For 
the sake of comparison it is also shown the values of Xmin 
and Q,n for the ACDM model. 

It is also instructive to evaluate the effect of adding 
the systematic uncertainties of the SNe analysis on our 
results. Considering only SNe, the addition of the 
systematic erros to the statistical erros lead us to flm = 
0.295^0036 foJ' the Visser's model. We also obtain a 
considerable lower value for the difference between the 
X^ of the two models Ax^ = 8.11. Now, taking into 
account the CMB and BAO measurements together with 
SNe, we obtain n„, = 0.290tJJ;^i^ and Ax^ = 10.26 (see 



FIG. 2: Best-fit for the distance modulus versus redshift for 
the Visser model (sohd gray fine) and the ACDM model 
(dashed line);_^The SNe data were taken from the Union2 
compilation [2^ . 



Table U. 

In the Fig. [1] and Fig. [2] we show the Hubble 
parameter and the distance modulus as functions of 
redshift considering the best-fit value of il,n for the SNe. 
For the sake of comparison, the standard ACDM model 
is also shown. Note that although the massive graviton 
model is disfavored, it seems to be able to reproduce 
very well the SNe la measurements, as can be seen in 
the Fig. [21 This shows the importance of the x^ test in 
distinguishing the two models. 



E. Effective equation of state 

The Fig. ([3]) shows the effective equation of state 

2(1 + z)dH 



Weff{z) 



-1 



iH dz 



(19) 



as a function of the redshift for the best-fit values 
above. The deceleration parameter, which is shown in 
the FigUl is related to Weff through q{z) — (Sweffiz) + 




FIG. 3: EfTective state parameter as a function of the redshift 
for the best-fit value obtained from SNe la. 




FIG. 4: Deceleration parameter as a function of the redshift 
for the best-fit value obtained from SNe la. Note that the 
acceleration phase is transient in the Visser model. 



V. CONCLUSIONS 

The theory of massive gravitons as considered in 
the Visser's approach has the advantage that the field 
equations ([3]) difi'ers from Einstein equations only in a 
subtle way, namely, by the introduction of the bimetric 
mass tensor M^^. Moreover the van Dam-Veltmann- 
Zakharov discontinuity (vDVZ) present in the Pauli- 
Fierz term can be circumvented in Visser's mode l by 
introducing a non-dynamical flat-background metric [30] . 

From the cosmological point of view, the meaning of 
the mass tensor, classically speaking, is a long range 
correction to the ordinary Friedmann equation. Such a 
correction mimics the effects of a dark energy component 
in such a way that additional fields are not necessary. 

In this context, we have shown that the cosmological 
model with massive gravitons could be a viable 
explanation to the dark energy problem. But, although 
the parameter 17^ is well constrained, the model 
is disfavored when compared to the ACDM model. 
Considering systematic errors, the difference between the 
X^in of the two models reduces considerably, but the the 
Visser model is still disfavored. 

Finally, the plots of the effective state parameter 
and of the deceleration parameter for the best fit value 
of f2m, show a very particular feature of the Visser's 
model, namely, the transient behavior of the accelerated 
phase of expansion. The Universe begins to accelerate 
approximately at the same redshift of the ACDM model, 
but for a very small redshift (2 ~ 4 x 10~^) we have a 
second transition and the Universe becomes to decelerate 
again. In spite of this, the behavior of the Hubble 
parameter H{z) is very similar in both models as can 
be seen in the FiglT] In this way, one would think that 
the transient acceleration phase is what make the Visser 
model less compatible with SNe data than the ACDM 
model. This is a problem which we will address in the 
future in order to find consistent modifications of Visser's 
approach. 



l)/2. In order to plot these curves we have included 
a component of radiation with the present value of 
the density parameter J7,. = 5 x 10~^. For the best- 
fit value found in our analysis, the Visser model goes 
through the last three phases of cosmological evolution, 
i.e., radiation-dominated {w = 1/3), matter-dominated 
{w = 0) and the late time acceleration phase {w < —1/3). 

Note that for low redshifts the Visser's model shows 
additionally a phase dominated by matter, indicating 
that for this model the late time acceleration of the 
Universe was a transient phase which has already 
finished. Moreover, for low redshifts, this behavior of 
the Visser's theory is in accordance with the fact that 
the luminosity distance of very low redshift SNe la can 
be fitted with CDM model only, i.e., at very low redshift 
the ACDM, CDM and Visser's model are degenerate for 
the cosmological observations. 
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